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The Fisher Matrix is the backbone of modern cosmological forecasting. We describe the 
Fisher4Cast software: a general-purpose, easy-to-use, Fisher Matrix framework. It is open source, 
rigorously designed and tested and includes a Graphical User Interface (GUI) with automated E^T^HjX 
file creation capability and point- and- click Fisher ellipse generation. Fisher4Cast was designed for 
ease of extension and, although written in Matlab, is easily portable to open-source alternatives such 
as Octave and Scilab. Here we use Fisher4Cast to present new 3-D and 4-D visualisations of the 
forecasting landscape and to investigate the effects of growth and curvature on future cosmological 
surveys. Early releases have been available at http://www.cosmology.org.za since May 2008 with 
750 downloads in the first year. Version 2.0 is made public with this paper and includes a Quick 
Start guide and the code used to produce the figures in this paper, in the hope that it will be useful 
to the cosmology and wider scientific communities. 

PACS numbers: 98.80.Es, 95.36.+x,02.70.Rr 



I. INTRODUCTION 



The need for rapid forecasts of constraints from proposed new surveys has played an important role in 
the approach of the cosmology community to the mysteries of dark energy. Faced with an almost total lack 
of understanding of the physics underlying dark energy, focus has shifted towards designing surveys that 
are optimal in some information theory sense. When hunting an entirely unknown animal, asking questions 
about its beak or tail may leave one barking up the wrong tree, so to speak. Instead a safe bet is simply to 
maximise the number of bits of information one gets about the animal. The Fisher Matrix formalism allows 
one to do this within a parameterised framework. 

Despite its limitations, the Fisher Matrix formalism has become the de facto standard in cosmology for 
comparing surveys and for forecasting parameter constraints^. It can also be easily adapted to different 
scenarios and has thus informed the design and funding of essentially all modern cosmological surveys. 

Nevertheless, while the Fisher Matrix formalism is relatively simple in principle, students can find it 
challenging at first, implementations are at times buggy and usually limited - almost always assuming a 
flat universe, for example - and there is unnecessary repetition of code, without a clear development path 
accessible to the entire community to build on. Fisher4Cast was initially developed in 2007 and 2008 to 
address these issues with the aim of providing a free, graphical framework that would make it easy for 
students to learn the formalism while also simultaneously providing a rigorous, robust and general code-base 
for researchers to build on and extend. Version 2.0, released with this paper, continues with the same aims 
as . 

We discuss the Fisher formalism in Section [Til and introduce the cosmology used in Fisher4Cast in Sec- 
tion [nil The effect of including cosmic curvature as a parameter is discussed in Section HV CI and the growth 
of structure as an observable in Section IIVBI The applications of Fisher4Cast are outlined in Section IIV Al 
while new features in Versions 1.2 and 2.0 are described in Section FlVDI The appendices give the general and 
explicit Fisher derivatives and the Quick Start guide for Fisher4Cast, as well as the samples of the Matlab 
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code used to produce the figures in this paper. 



II. FORECASTING AND THE FISHER MATRIX 



Forecasting of survey constraints for a proposed survey can be achieved either through full Monte Carlo 
simulations of the survey which are typically time consuming and computationally intensive, or by using the 
much simpler and quicker Fisher Matrix technique, which wc now discuss in detail and which forms the basis 
of Fisher4Cast. The Fisher Matrix translates errors on observed quantities measured directly in the survey 
into constraints on parameters of interest in the underlying model. Put more directly, it is the elegant way 
of doing propagation of errors in the case of multiple, correlated, measurements and many parameters [4]. 




As an example, consider an arbitrary function y = f[z, 0) of some parameter 9 and an independent variable 
z. Assuming a perfect measurement of z, the error 66, for a given measured Sy is, by simple calculus, 
SO = {^)~^Sy, or equivalently {36)^'^ = (§9)^(^2/)"^- This is perhaps the simplest example of a Fisher 
Matrix; with a single element, cf. Eq ([3]). Here 6 represents the parameter we want to measure and / the 
observable quantity (e.g. H{z)). 

In more generality, the Fisher Matrix formalism predicts the constraints on a vector of parameters 6 ~ 
(01, 02, ■■■10 A, •■) - such as wq and Wa - resulting from measurements of one or more observablcs X" = X"(0, z) 
(such as H{6,z) or (i^(0, z)), each at a range of redshifts, z = {zi, Z2, ■■■, Zi, ..) e.g. in a BAO survey one 
might measure H{9,z) and dA{0,z) at a single rcdshift, while a Type la supernova (SNIa) survey may 
measure dL{0,z) at hundreds of redshifts. There arc therefore three indices in general to keep track of, 
(A, a, i) corresponding to parameter, observable and rcdshift. The number of observablcs is arbitrary and 
combining results from independent observablcs is essentially trivial, hence we will often suppress observable 
index, X" = X for simplicity. Boldface indicates the entire vector, either of parameters 6, observablcs X or 
redshifts z. 

The Fisher Matrix estimates not only the individual errors on the parameters, 9, evaluated at a given 
input /base/fiducial model 6 — 6*, but also the correlations between them, leading to the characteristic 
Fisher error ellipsoids (ellipses if one considers only pairs of 9a)- To make this clear, consider the likelihood, 
^ = P{d\9), for a general survey, which gives the conditional probability of observing the data d = 
(di, (i2, rfj, ■•) assuming the cosmological model 6 is correct. We can expand the likelihood around the 
fiducial model: 



where Da = 99 a represents the partial derivatives with respect to the parameter 9 a. The first term in the 
expansion is a constant depending on the fiducial model. The fiducial model is expected (after averaging over 
many data realisations) to be the point of maximum likelihood, hence the first derivative of the likelihood 
vanishes. The third term is the curvature matrix/Hessian of the likelihood, and is the term used in the 
Fisher Matrix which is formally defined as the expectation value of the derivatives of the log of the likelihood 
with respect to the parameters 6, or 



Here the angle brackets indicate the expectation value which, for an arbitrary function g{;x), is defined to 
be {g{X)) = g{x)fx{x)dx with fx{x) the probability distribution function of the random variable x, 
which here is the noise on the data, assumed to be Gaussian of mean zero. 

The likelihood for a given observable X is expressed in terms of the theoretical value of the observable 
Xi evaluated at the redshifts Zi and data for that specific observable as .if oc exp(— A-'^C^^ A/2) where 
A = X — d, generalising the usual chi-squared statistic relating the theory by allowing for a general data 
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FIG. 1: The link between model curves, data and error ellipses. Parameter values inside the ellipse are 
a good fit to the data, such as the blue dot-dashed curve (left) corresponding to the open blue circle (right). The 
horizontal line dA — dA.ACDM = corresponds to the fiducial model on the right {wq — l,Wa = 0)- A parameter value 
outside the ellipse is not a good fit to the data however, as illustrated by the orange dashed line (left) and orange 
cross (right). The ellipse and simulated error bars in the two panels both correspond to a nominal survey measuring 
the angular diameter distance, dA{z), with 6% error in ten bins between 0.1 < 2 < 3, with ACDM assumed to be the 
correct model in both cases. 



covariancc matrix C. Substituting the above expression into Eq. ([2]) converts the equation from derivatives 
of the likelihood itself into a sum over derivatives of the observable X with respect to the parameter 9a' 

i * 

where dC/dOx is the derivative of the data covariance matrix with respect to the parameter Ox which is 
assumed to vanish in the second equality implying that the data errors are independent of cosmological 
parameters. This is often the case, e.g. the errors on measurements of Type la supernova (SNIa) flux are 
independent of the dark energy parameters Wi^^Wa to good accuracy. The second equality also requires that 
the data are uncorrelated, in which case C is diagonal with entries erf, with the Ui the 1 — cr error on the 
i-th data point. 

In the case where we have multiple independent measurements of different observables X" (e.g. H{z) and 
dA{z)), the total Fisher Matrix is just the sum of the individual Fisher matrices indexed by a. Similarly, if we 
have independent prior information, this is encoded in a prior matrix between the cosmological parameters. 
In this paper we will refer to the prior on a single parameter 6x as Prior(0A) = (A^a)"'^, where A^a is 
the uncertainty on the parameter as measured from prior surveys, see e.g. Figure (jlip . In the case where 
the different measurements are not independent, they must be combined with the suitable data covariance 
matrix. The inverse of the Fisher Matrix, i^AB' provides an estimate of the error covariance matrix for the 
parameters 0a, as we now expand upon. 

In general one often considers a p + 71-dimensional likelihood which includes not only the p parameters 
of interest but an additional n nuisance parameters that form a natural part of the problem but are of no 
direct interest (such as Hq in studies of dark energy dynamics). As a result we are usually only interested in 
the likelihood as a function of the p key parameters, which is obtained by marginalising over the nuisance 
parameters, viz. 



p) = / ^{9i, ..,9p, ..,9p+n)d6p^i...d9p+n- (4) 



For a likelihood of arbitrary shape marginalisation must be performed numerically and is well-suited to 
Markov Chain Monte Carlo (MCMC) methods but in the special case where it is a multivariate Gaussian, 
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FIG. 2: Varying wo and Wa of the fiducial model — the generated ellipses for a measurements of the Hubble 
parameter H{z) and angular diameter distance dA{z) survey characterised in Table U The values for the coefflcients 
in the CPL parameterisation for w{z), wo and Wa, have been varied on a grid over —1.3 < wo < —0.6, —0.7 < Wa < 1- 
As wojWa change, not only does the ellipse centre shift, but the size of the ellipse changes, as well as the slope of the 
degeneracy direction between the two parameters. Since the Dark Energy Task Force Figure of Merit is linked to the 
inverse of the area of the ellipse the value of the FoM increases for with increasing wo and Wa- 

the marginalisation can be performed analytically and simply in terms of the Fisher Matrix. Let us write 
the full Fisher Matrix in terms of sub-matrices, as 



where is the p x p sub-matrix corresponding to the parameters of interest, O is an p x n matrix describing 
the correlation between the nuisance parameters and the parameters of interest'^ and Af is the n x n matrix 
representing the nuisance parameters we wish to marginalise over. The marginalised Fisher Matrix for the 
parameters of interest is then given by 0] : 



where the first term is the matrix of the reduced parameter space of interest, the second term encodes the 
effects of the marginalisation over the other nuisance parameters and ^ represents matrix transpose. 

As we mentioned before, the inverse of the Fisher Matrix provides an estimate of the parameter covariance 
matrix. For an unbiased estimator (that is one whose expected value of is equal to the fiducial model 
6* assumed to be correct), and in the case where one does not marginalise over any other parameters (i.e. 
we consider all other parameters perfectly known), the expected error on any parameter satisfies the 
Cramer-Rao bound^ 

A0A > , (7) 

while in the more realistic case that one wants to marginalise over all the other parameters in the problem, 
the bound becomes 




(5) 



F = e- OAf-'^o'^ 



(6) 



A0A > V(F-i)aa 



(8) 



^ If O = then the nuisance parameters have no impact on 0. 

* For a proof see e.g. , http: / / en.wikipedia.org/ wiki/CRLB, or e.g. p. 426 in Q 
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i.e. one first inverts the Fisher Matrix, then takes the AA component of the resulting matrix. One can 
show that the latter is always greater than or equal to the former, i.e. marginalisation cannot decrease the 
error on a parameter, and only has no effect if all other parameters are completely uncorrelated from the 
parameter of interest. Note that in the case where the likelihood is exactly Gaussian in the parameters, the 
Cramer-Rao bound becomes an equality and not just a lower bound. 

Since Fisher Matrix analysis assumes the likelihood is a multivariate Gaussian (an approximation that 
can be made arbitrarily good by considering better and better surveys), contours of constant probability are 
ellipsoids within the Fisher formalism. These ellipsoids (ellipses for two parameters) are given by solving the 
equation 

A6»^FA6/ = /3 (9) 

where = 9 — 6* is the parameter vector around the fiducial model, 6* , and /3 is a constant determined 
by the desired confidence level and the number of parameters. For two parameters, the 1 and 2 — cr contour 
levels correspond to /3 = 2.31 and 6.17 respectively 0, The Fisher4Cast GUI allows plotting of both 1- 
and 2-dimensional contours and hence always marginalises the full 5-dimensional Fisher Matrix to achieve 
this. Marginalisation over some or all of the other parameters can be effectively switched off by making the 
corresponding diagonal elements of the prior matrix very large. The Fisher Matrix and the corresponding 
ellipses provide the Gaussian estimate for how well the parameters of the model will be constrained by a 
given experiment assuming the true model is that at which the Fisher Matrix was evaluated (e.g. ACDM). 
This is illustrated in Figure H]), which shows the 1 — cr error ellipse around the fiducial ACDM model with the 
coefficients in the Chevallier-Polarski-Linder (GPL) (see Eq. (fTTj) ) parameterisation {wq, Wa) = (—1, 0), 

for a survey consisting of measurements of the angular diameter distance between redshifts of 0.1 and 3. 
Values of wojWa inside this ellipse will have expected likelihoods that differ from the fiducial model by less 
than 1 — (7. The Fisher Matrix allows us to estimate which sets of parameter values wc will be able to rule 
out at a given significance level if the fiducial cosmological model is correct. Wc will see in Section IIV Al 
(see Figure ([5])) that changing the assumed fiducial/base model has a big effect on the ellipses for the same 
survey. 

While the ellipses provide significant insight they do not allow immediate comparison between different 
surveys, a feature required if one wants to optimise or compare surveys head-to-head [sl. [ol. [lOl [III. [T^. [isl 
[Tsj . One common way to perform such comparison is to formulate a Figure of Merit which ascribes a single 
real number to each survey, the simplest of which is to use the volume of the ellipsoid or of the marginalised 
ellipse. The volume of the n-dimcnsional error ellipsoid (corresponding to an n-dimcnsional Fisher Matrix) 
is: 

/3 



where P is defined in Eq ([5]), Vsn = 7r"/^/r(^ + 1) is the volume of the n-dimcnsional unit sphere and T{u) 
is the Gamma function. For the interesting case n = 2, Vga = tt is, of course, the area of the unit circle. 
Note that it is common in the literature to ignore the Vs" and /3 factors and to incorrectly refer to the 
determinant factor alone as the area or volume of the ellipse/ellipsoid. 

In fact, since the Fisher Matrix is a metric, the square root of the determinant is a natural volume element 
providing the Jacobian for the action of the linear mapping induced by the Fisher Matrix, i.e. one should 
think of the Fisher Matrix as inducing a linear mapping rather than 'being' an ellipse itself. 

Fisher4Cast includes the standard FoMs as well as some new ones available through the GUI and command 
line. Those using the volume are based on Eq PU)) with n = 2. Although some of the Figures of Merit are 
only defined for the error ellipse in the wq — Wa plane, where wq, Wa are the coefficients in the GPL 0, 01 
parameterisation of the equation of state of dark energy (see for e.g. [i3l), the FoMs in Fisher4Gast are 
calculated by the code for any pair of cosmological parameters being considered rather than the full 5-D 
matrix. We briefly outline the FoMs used in Fisher4Gast: 

• DETF 

This Figure of Merit in the Report of the Dark Energy Task Force [l3| is defined to be the reciprocal of 
the area of the 2 — a error ellipse in the wq — Wa plane of the GPL dark energy parameterisation 0]. 
This is, via Eq (fTO|) . equal to det(F^/^)/(7r-\/6.17). Unfortunately the DETF report does not appear 
to use this definition, and instead quotes det(F^/^), which is the inverse of the 1 — a ellipse in units of 
the area of the unit circle. Because of the benefits of the geometric interpretation Fishcr4Gast returns 
the true inverse area of the 2 — ct ellipse. To convert from one DETF FoM to the other, one should 
muhiply the Fisher4Gast DETF output by nVKlT ~ 7.8. 
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• Area]^_g. 

This Figure of Merit is the reciprocal of the 1 — cr error elhpse area in the parameter plane currently 
plotted, i.e. det(Fi/2)/(7rV2:3T) 

• Areai_o- 

Simply the inverse of the previous FoM. 

• TrC 

This FoM is defined as the trace of the covariance matrix of the data, C = F^-*-, estimated as the inverse 
of the marginalised Fisher Matrix. This FoM is simply the sum of the squares of the marginalised errors 
on each parameter. 

This FoM is defined as the sum of the squares of the entries of the whole covariance matrix, C = F"-*-. 
Unlike the previous definition this FoM is sensitive to the off-diagonal components of the covariance 
matrix as well as the diagonal components. 

Designing a FoM which is both robust and easy to use, while not favouring a particular model space or 
paradigm has been the topic of recent work [l6|. Other suggestions have been proposed, such as using 
principal component analysis , as the basis of an effective FoM to compare survey efficiency. 

Finally we note that in addition to the Fisher Matrix, one can include the non-Gaussian terms in the 
expansion of the log likelihood (Eq. ([T])). Study of these flex corrections is left to future work. Their inclu- 
sion in Fisher4Cast is currently implemented and will become public with the next version of Fisher4Cast, 
providing a test of the Fisher Matrix approximation. 



III. THE COSMOLOGY OF HUBBLE, DISTANCE AND GROWTH 

Although Fisher4Cast is a completely general Fisher Matrix framework at the command-line level, the 
GUI is coded as a cosmology interface, since this is its primary application. In the context of modern 
cosmological surveys, the primary observables are the expansion rate of the Universe, measured through the 
Hubble rate H{z), cosmological distances such as the angular diameter distance, ^^(z), and the growing 
mode of dark matter density perturbations, 5(x, z) oc G{z). H{z) and dA{z) are provided by BAG surveys 
while growth can be measured using lensing or number count surveys and potentially also BAG if the bias 
is measured independently (e.g. through redshift distortions - see [isl [l9|). 

The Fisher4Cast GUI uses the observables if, (Ia and G in a general Friedmann-Lemaitrc-Robcrtson- 
Walker (FLRW) universe. The cosmic parameters assumed for the GUI are (iJo, i^m, ^^fc7 "Wo, li'a), where i/o 
is the value of the Hubble constant in kms~^Mpc~^, fim is the energy density of matter today in units of 
the critical density, Q,k is the curvature energy density (JIde = 1 — — ""^O) '"'a s-^^ the coefficients 

in the GPL expansion of the dark energy equation of state [1, 0] : 

w{z) = Wo + Wa-^— = Wo + Wa ( 1 - — ) , (11) 

1 + z \ ao J 



where ao = c/(i7o v|f2fe|) is the curvature radius of the cosmos. Gthcr dark energy expansions can be easily 
accommodated in Fisher4Gast by changing the appropriate input functions. The expansion history of a 
FLRW universe is described by the Hubble parameter: 

H\z) = H^E\z) = Hi + zf + 17fe(l + zf + (1 - an - a)/(^, wo, wj) , (12) 

with the evolution of the dark energy density, jOde(z) oc f{z) determined by 

nz) = e.^(,p-±^dz). (13) 



^ Unless explicitly indicated elsewhere, references in this paper to Hi will mean the current value of the density parameter and 
not its value as a function of time. 
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FIG. 3: Fisher Derivatives for the Hubble parameter (left panel), angular diameter distance (middle panel) and 
growth function (right panel) are shown for the parameters considered in our cosmological example: Hq (solid dark 
red line), In Sim (dot-dashed red line), 51^ (dashed dark orange line), wo (dotted orange line) and Wa (dot-dashed 
peach line). The full set of analytical derivatives of the H{z) and dA{z) are found in AppendixlB] In the case of the 
growth function G{z) the derivatives of the solution to Eq. (|21|l are taken numerically. 



For the CPL parameterisation, Eq. pT|) . f{z) is given by 



f{z) ^ + 



exp 



1 



(14) 



The angular-diameter distance, dA{z) relates the angular size of an object to its known length, providing 
a measure of the distance to the object, and is given by: 

1 



dA{z) = 



where 



1 + zHo 



X{z) = 



sinh ( y 



dz' 



^kXiz] 



(15) 



(16) 



/o Eiz')' 

and E{z) is as defined in Eq. ([T^ . These forms are valid for all values of fJfc via continuity and the 
trigonometric identity sinh(ia::) = i sin(a;). The often-used equation for the angular diameter distance contains 
three equations, depending on the sign and magnitude of flk, however this is redundant, at least conceptually. 
In numerical analysis we use the Taylor series expansion for very small flk, Eq. (jBlOp . 



ddAiz) 



an. 



1 



Hol + z 



where 



X{z,0)^ Xiz)\,,^_ 
x{z)) assuming flatness. 



(17) 



(18) 



are the functions (for example E{z) 

Finally we discuss the gove rning e quat ion for the growth of structure, a potentially powerful probe of dark 
energy [l3, [10, [U, [13, [UHli, [H, (li, M [H, HI. In general one needs to solve the differential equation for 
the perturbations in the matter density 6 (assuming the pressure and pressure perturbations of the matter 
are zero - p = Sp = 0) [13, [Ml, [13 : 



(19) 



We discuss this equation in the context of curved universes with dynamical dark energy in Section IIV Bl 
Fisher4Cast takes as input constraints on the growth G{z) which provides the temporal evolution of density 
perturbations, i.e. 6{x.,z) oa G{z). 

Error ellipses for all the cosmic parameters in our example are shown in Figure ([3]), where they have 
been computed around the flat-ACDM fiducial model - (ifo, ^^m, ^fc, wo, Wa) = (70,0.3,0,-1,0) - and the 
full set of analytical derivatives for H{z), dAiz) are given in the Appendix[Bl In the case of the derivative of 
the angular diameter distance in terms of the curvature parameter ddA/dilk, which will be zero in the flat- 
ACDM model, a Taylor expansion of the derivative (Eq. [B9|) around il^ — is also provided in Section IB] 
The derivatives for growth are computed numerically since no general analytical solution for the growth 
exists; see Section llVBI 
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FIG. 4: Marginalised Fisher ellipses for all parameters. Fisher error ellipses for the full range of cosmo- 
logical parameters, namely Ho, ilm,^k, wo, Wa, (for a survey characterised by Table U with the addition of growth 
measurements at the same redshifts and with the same precision as the Hubble parameter) are produced by calling 
Fisher4Cast in a simple loop over all parameters. The dark and light filled ellipses indicate the 1— and 2 — a contours 
respectively marginalised over all other parameters. The diagonal panels show the fully marginalised one-dimensional 
likelihood for each parameter; 1 and 2 — a limits shown as solid and dashed vertical lines respectively. The code used 
to produce these plots is included in the latest release of Fisher4Cast. 



IV. FISHER4CAST AND ITS APPLICATIONS 

Fisher4Cast is written in Matlab^ following an object-oriented model and using standard software engi- 
neering standards for implementation and testing. The code is not specific to cosmology; Fisher matrices 
can be generated given any X(0). The Fisher derivatives are computed analytically (if they are 

known) or numerically, as in the case of growth, allowing the code to handle complex cases without ana- 
lytical formulae for X. The code suite includes a Graphical User Interface (GUI), which is specific to the 
cosmological example we discussed above in Section Hill Fisher4Cast includes an automated M^rpC report 



See |http:// www. mathworks . com I 
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generating feature described in detail in the Users' Manual'' [33|, which automatically creates a summary of 
relevant data, matrices and figures. 

Fisher4Cast facilitates novel research and education in two different ways. Apart from being well-tested 
against existing Fisher Matrix results, it is general and modular. Because of this modular nature, a natural 
application of Fishcr4Cast is to visualisation. The code can easily be called repeatedly in large loops, enabling 
one to study large-scale properties of the Fisher Matrix for a wide range of surveys and cosmologies. We give 
examples of such studies in the following subsection. Secondly, Fishcr4Cast is coded for a general FLRW 
universe (the curved case has rarely been studied in the literature). We study the issues of curvature and 
growth in dark energy Fisher analysis in the last two subsections. 




FIG. 5: Figure of Merit plane — the Dark Energy Task Force Figure of Merit (FoM) for a survey consisting of one 
bin each of the Hubble parameter and angular diameter distance, with the fractional errors an /H = adj^/dA = 0.1. 
The redshifts of the H,dA measurements are varied separately from z = 0.1 to z = 5, and the resulting FoM for 
each survey configuration is plotted as a 3-dimensionaI landscape (left panel), or a flat 2-dimensional plane (right 
panel). The colourmap in both panels reflects the value of the FoM, from low values of FoM ~ O(blue) to higher 
values (FoM ~ 0.08) . 



A. Visualisations 

One of the key design principles behind Fisher4Cast is ease of use - the GUI was specifically created to 
provide users without an in-depth knowledge of Fisher Matrix theory access to the power of the formalism. 
Fisher4Cast can also be called from the command-line, as shown in the code examples in Appendix[Cl In this 
subsection we investigate applications of Fisher4Cast to probe and visualise the Fisher Matrix. Figure ([2|) 
illustrates the resulting ellipses when Fisher4Cast is called in a simple two-dimensional loop that varies 
the values of wq — Wa in the cosmological model assumed to be true. The strong dependence of both the 
orientation and size of the resulting ellipse on the assumed model is clearly evident, as the ellipses rotate 
and shrink for larger values of wq and Wa- 

Figure (O illustrates the complementary case when the assumed cosmological model is kept fixed, but 
where instead the survey parameters are varied. In all cases the survey consisted of one measurement of dA{z) 
and H{z) in a single redshift bin (with fractional errors of 10% on cither observable), while the redshifts of 
these two bins were varied independently over the range 0.1 < z < 5. The DETF Figure of Merit (FoM) 
of the survey (cx: 1 /Area of the wq — Wa ellipse) is plotted as a FoM "landscape" , where the colourmap is 
related to the value of the FoM - the higher values of the FoM are depicted in red, and the lower values in 



The Users' Manual is also bundled with the Fisher4Cast code software available at [s^l . 
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blue. The FoM ranges between 10~^ and ~ 0.08. The right-hand panel of Figure ^ shows this landscape 
projected into two dimensions, with the same FoM colourmap as in the three dimensional case. 

Two interesting features are immediately apparent. First, the peak in the FoM landscape, or "FoM 
hotspot" occurs for a survey with measurements of H(z) at z ~ 1.5 and dA{z) at z = 0.6, while a ridge 
of moderately good values of the FoM (relative to the average) emerges along the line corresponding to 
H{z — 0.2) for measurements of cIa and redshifts larger than 1.6. Secondly these two regions of higher 
Figure of Merit are separated by a "cold valley" of lower FoM values. This is of particular interest given 
the optimisation of current and future Baryon Acoustic Oscillation (BAO) surveys [i,i,[i0„iii,^- Future 
BAO surveys will measure the radial and tangential oscillation scale (see [3a] for a recent review of BAO in 
cosmology), producing measurements of cIa and H at the same redshifts, i.e. along the diagonal line in this 
plot. For z roughly between 0.5 and 3, this line intersects regions of relatively high FoM. This landscape is 
explained by the geometric interplay between the size and orientation of the degeneracy directions between 
the H{z) and dAiz) ellipses, as is illustrated in detail in Appendix [XI Finally we extend the H — cIa landscape 




FIG. 6: 4D Fisher Hypersurface Visualisations — DETF FoM for the combination of a single measurement each 
oi H, dA and G (where the growth is normalised to unity at z = 0). The redshift of each measurement is allowed to 
vary, generalising Figure ((5)1 by including growth. Slices through this hypersurface (left-hand panel) at zh = 1.67, 
Zdji = 0.67 and zg = 2.50 show the relationship between the redshifts at which the measurement of the observables 
are made and the FoM. The colourmap goes from a FoM of 9.5 x 10^'' (dark blue) to 0.14 (dark red). Adding a 
high-redshift measurement of the growth function tightens the constraints on the dark energy parameters wo and Wa , 
shown by the red ridge of high FoM values. This hypersurface can be illustrated in a complementary way: surfaces 
of constant FoM are shown in the right-hand panel - ranging from 0.03 (transparent light blue outer surface) to 
0.15 (dark red opaque centre surface). Comparing the left and right panels, one notes that the intersection "hotspot" 
region in the left-hand panel at {zh = 0.4, Zd^ ~ 0.67, zg ~ 5) is contained within in the red iso-surface with the 
high value of the FoM. This iso-surface extends right down to redshift zg ^ 1, re- iterating the improvement on dark 
energy constraints when including growth. 

to include a loop over the redshift of a single growth measurement in the range 0.1 < z < 5. In all cases we 
normalise the growth today; G{z = 0) = 1. The full four dimensional surface cannot be plotted in general, 
we show slices through the hypersurface in Figure ([6]), illustrating how a measurement of the growth at 
high-redshift measurement leads to much higher values of the FoM overall, and opens up interesting new 
"hotspots". These are only a few of the potential visualisation applications for Fisher4Cast. 

B. Dark Energy Constraints from Growth of Structure 

The growth function G(z), defined as the solution to Eq. ([T9|) is sensitive to dark energy; the Hubble pa- 
rameter acts as a friction term in the differential equation, increasing or suppressing the growth of structure. 
While in general there is no analytical solution to Eq. (|19p . under the assumption of a flat universe and a 
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cosmological constant (or pure curvature) the growing mode satisfies the following integral form [3y, [sj : 
where the 5/2 coefficient is chosen to ensure that G{z) 1/(1 + z) during matter domination. This 
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FIG. 7: Comparing analytical derivatives to the full solution of the ordinary differential equation. The 

left-hand panel in this plot shows the derivatives of the growth function G{z), where analytical derivatives are taken 
of Eq. (|20[l dashed red lines, or where derivatives have been taken numerically of the solution to Eq. (|21|l solid brown 
lines. The parameters considered are wq (light red for the analytical derivatives and dark brown for the numerical 
derivatives) and Wa (dark red and light brown for analytical and numerical derivatives respectively). The right-hand 
panel shows the ellipses that result from a survey of the growth measured in 20 redshift bins from 0.1 to 2, with 
10% error on the growth function (normalised to unity at z = 0), and priors of lO** on the Hubble parameter, matter 
and curvature densities respectively. The ellipse that results from the analytical derivatives red dashed line suggests 
(incorrectly) much tighter constraints on the dark energy parameters, whereas the ellipse from the numerical solution 
(solid brown line) is much more degenerate in the dark energy parameters. 



Parameter 


Value 


Redshifts of 


H : z = [0.3,0.6,0.8,1.0,1.2,3] 


measurement 


dA-. z^ [0.3, 0.6, 0.8, 1.0, 1.2, 3, 1000] 


Percentage error 


Value [%] 


Hiz) 


an/H =[5.80, 5.19, 3.59, 2.84, 2.53, 1.48] 


dA{z) 


(Tdyi/dA = [5.19, 4.30, 3.22, 2.3, 2.03, 1.19, 0.22] 


Cosmological model 


Value 




(70kms-^Mpc-\ 0.3, 0, -1, 0) 


Priors on model 


(lO"*, 10^, 10^, 0,0) 



TABLE I: Survey data from the Seo & Eisenstein survey configuration jSSi ] — used in Figures ([8]), (fTOjl . ^ 

and Q. In some cases measurements of the growth function were added, taken at the same redshifts as the Hubble 
parameter; in others the prior information on various parameters was changed. See the captions of the relevant figures 
for the specific details. 



expression should not be used however, to compute the Fisher derivatives dG/dVLk^dG/dwo or dC/dwa 
since all of these derivatives violate the validity of the equation. Instead, the growth derivatives should be 
computed numerically from the solution of the full differential equation for 5{x). Rewriting the Raychaudhuri 
equation in terms of the Friedmann equation and the curvature density allows one to find an equation 
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FIG. 8: Precision measurements of the Growth tighten dark energy constraints — we show how changing 
the error on the growth relative to the Hubble parameter, 5G/G = aSH/H, from a = 0.1 (innermost ellipses) to 
« = 10 (outermost ellipses) affects the power of the survey. We consider the survey summarised in Table U] but 
with the addition of growth function measurements with 10% error evaluated at the same redshifts as the Hubble 
parameter, for simplicity. The left-hand panel shows the Fisher ellipses in the wo — Wa plane, while the right-hand 
panel shows the DETF FoM as a function of a. As the error on growth is decreased towards very small values 
(crG/G 0.1 — an/H), the FoM increases as expected. For large values of a, however, the FoM flattens out as there 
is essentially no pertinent information from the growth function. 



explicitly showing the curvature and dynamical dark energy contributions to the friction term: 



where the new independent variable is a; = a/ao ~ 1/(1 + z), oo is the radius of curvature and Qk{x) ~ 
— fc/(aQa;^i7(x)^); f2DE(a^) = Pde{x) / Pcrit{x) are the fractions of the critical density in curvature and dark 
energy respectively. Alternatively, this can be written as a differential equation in terms of ln(x): 



which is the equation actually solved in Fishcr4Cast since it is typically more stable numerically. Appropriate 
initial conditions for this differential equation arc set deep in the matter dominated era: 
G{zi) = l,dG/d\nx{zi) = G{zi) for Zi > 100.^ Note that as a result, the growth solutions will be unreliable 
if w{z —^ oo) = Wo + Wa > (or even if it is close to zero from below) since then there will be significant 
or even dominant early dark energy. Fisher4Cast allows the user to choose the redshift where the growth is 
normalised to unity. The Fisher derivatives all satisfy dG/dOi = at the normalisation redshift. 

While the growth functions Eq. (|20p and Eq. agree for ACDM, analytical Fisher derivatives taken 
of Eq. (|20[1 will not agree with the numerical derivatives of Eq. (|22p . and will thus produce very different 
error ellipses when used incorrectly. The left-hand panel of Figure ^ shows the wq, Wa Fisher derivatives of 
both solutions for ACDM. The derivatives are overestimated for both parameters, leading to (spurious) tight 
constraints on the parameters in the Fisher ellipse, as is illustrated in the right-hand panel of Figure ([7]), 
for a survey characterised by 20 measurements (at 10% accuracy) of the growth function between z = 0.1 
and z = 2. This illustrates the danger of using derivatives based on the analytical form of G, even when 
evaluated at the model for which the analytical form is itself valid. 

As the precision of growth measurements increases, it is natural to ask how this impacts constraints on 
dark energy. This is easily investigated in Fisher4Cast. Consider Figure ([5]), which illustrates constraints 
on the CPL parameters around ACDM for a survey characterised by Table [H with the addition of growth 



One can compare with the growth code at |http: / / gyudon.as.utexas.edu/ ~komatsu] 




(21) 




(22) 
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measurements at the same redshifts as the Hubble parameter. As the error on growth relative to the error 
on H decreases, the Figure of Merit of the survey increases dramatically, highlightin g th e merit in ma king 

t Imprecision measurements of the growth function, as is the focus of recent interest |2ll . [H, [2^ [H, [l^, [281 
[3i,|40,[4l|,[42. 
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FIG. 9: Normalisation of the Growth function changes constraints — the wq — Wa error ellipse for mea- 
surements of the growth combined with one measurement of H at low redshift. The normalisation redshift of the 
growth function is varied from 2 = (dot-dashed line) to z — 4 (solid central black line) to illustrate the change 
in the degeneracy direction of the dark energy parameters with growth normalisation. Note that the constraints do 
not scale monotonically with redshift: the ellipse is most degenerate for normalisation of the growth at a redshift of 

2=1. 

As a final remark, we show how the redshift at which one normalises the growth function (normalisation 
redshifts z = 0, 1,2,4 are shown) influences dark energy constraints in Figure The constraints on the 
dark energy parameters are plotted for 5 measurements of the growth function at z = (0.3,0.6,0.8,1.2,3) 
each with a 10% accuracy, and one measurement of the Hubble parameter at z = 0.3 with <jh /H = 5.19%. 
An interesting point is that the resulting dark energy constraints do not scale monotonically with the 
normalisation redshift. 



C. The Effect of Cosmic Curvature On Dark Energy Constraints 

The degeneracy between curvature and dark energy has been well-studied even for perfect measurements 
of any single observable such as H{z) [i^, [3, |4^, |46|, [l^l, implying that marginalising over the curvature 
is important when performing parameter estimation and forecasting of constraints on dark energy. The 
degree to which curvature affects dark energy constraints is shown here as a simple example of FisherdCast. 
Figure (jlOp shows Fisher error ellipses for the dark energy parameters wo,Wa, after marginalising over 
curvature, as the prior information on curvature is changed from Prior(rife) = 10 (weak) to Prior(r2fc) = 
10^ (almost perfect) for the observables H, dA and G considered separately and in combination. While 
uncertainty in the curvature of the universe (represented by a small prior value on Qk) degrades all ellipses, 
this is much less pronounced when the observables are considered in combination, showing the importance 
of combining multiple probes of dark energy. For each of the parameters. Figure (|11[) shows that the 
constraints are eroded and the ellipse increases in size with the decrease in the prior - which expresses our 
confidence in the flatness of the universe. While the ellipses from the single observables such as H, dA show 
a much greater increase in size with decreasing curvature prior, using a combination of parameters is more 
robust - since combining multiple probes helps break the curvature-dark energy degeneracy. The change in 
the size of the ellipse directly relates to a change in the Figure of Merit (FoM), or power, of the particular 
combination of observables. Various FoMs are described in Section [Hi In the case of the FoM of the Dark 
Energy Task Force (DETF) [l3|, the area of the ellipse and the DETF FoM are inversely proportional. The 
right-hand panel of Figure (jlip shows this change in the FoM with changing curvature prior. Interestingly, 
the DETF FoM flattens out for both very large and very small values of the prior. This can be explained 
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FIG. 10: Curvature marginalisation and dark energy constraints - Fisher4Cast ellipses on the wq — Wa plane 
where the prior on the curvature density is changed from very strong, with Prior(nfe) = 10^ (left-hand panel), to 
very weak with Prior(f2fe = 10) (right-hand panel). As the prior is decreased (i.e. we are less confident in the value 
of ilk), the ellipse in the wo — Wa plane enlarges, in turn reducing the Dark Energy Task Force Figure of Merit (for 
the inner blue ellipse) from 3.298 to 1.890. The fiducial model is flat and the ellipses correspond to measurement 
of (moving from outwards in): G{z) (light brown outer band), H{z) (dark brown filled ellipse), dA{z) (orange filled 
ellipse) and the combination of all three (blue filled inner ellipse, outlined in black). The surveys details are given in 
Table HI with the addition of 10% measurements of the growth at the same redshifts as those of H{z). The prior on 
the matter density is kept fixed at 100. 



by the fact that for very large prior values, marginalising over the curvature has no effect, since our error on 
the curvature is minute, whereas for very small values of Prior(rii.) ~ the ellipse is completely degenerate 
with the curvature and hence reaches a maximum size. 



D. New Features in Fisher4Cast 

Fisher4Cast made its public debut (Version 1.1) in May 2008, with significant updates and revisions 
every six months since then. This work coincides with the release of Version 2.0 of the code. The main new 
features of versions 1.2 and 2.0 are listed below; a more comprehensive description of the changes in the new 
code is contained in the Readme.txt file, while the new features are described in the Fisher4Cast Users' 
Manual [s^ , provided with the distribution of the code suite. 

• Report Generating Features 

Code is included in the current version for automated generation of both I^TfTpC and text reports 
containing information such as the input survey, fiducial cosmology, output Fisher matrices and Fisher 
ellipse figure. 

• Fitting formulae for Baryon Acoustic Oscillation (BAO) Surveys 

Two extensions are provided to calculate the errors on the BAO oscillation scale according to the 
prescriptions of Blake et al. [1^ and Seo and Eisenstein [i^ (see the files in the directories called 
EXT_FF_Blake_etal2005 and EXT_FF_SeoEisenstein2007). These are not available via the GUI 
currently, and must be run directly from the command line. However, these modules combined with the 
rest of Fisher4Cast provide the capability of going directly from BAO survey specifications (volume, 
area, number density etc..) to dark energy constraints. 

• "Point- and- click" ellipse plotting 

A new feature allows the user to click in the figure and have the ellipse automatically generated at that 
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FIG. 11: Curvature marginalisation degrades dark energy constraints - using tiie same input survey data 
as Figure (|10|l . tiie error ellipses in the dark energy parameters for prior values on fife from l(f (inner curves) to 10 
(outer curves) are shown. Ellipses are produced for measurements of (clockwise from top left) the Hubble parameter 
(dark brown curves), the angular diameter distance (orange), the growth function (light brown), a combination 
the Hubble parameter and the angular diameter distance (green) and a combination of all three observables (blue). 
The bottom right panel shows the DETF Figure of Merit [10] (the inverse of the area of the ellipse in the uiq — Wa 
plane) for these ellipses as a function of the prior on Q.k- The colours correspond to the ellipses: green (dashed) is H 
and cLa combined while orange (dot-dashed) corresponds to (Ia alone. In all cases the shape of the curves is roughly 
the same, with no change in the FoM beyond a certain value of the prior, both in the case of very large and very small 
priors. Interestingly the curves for the Hubble parameter (brown) and angular diameter distance (orange) cross - a 
survey consisting only of measurements of H may seem to yield tighter (or weaker) constraints on dark energy than 
a survey only of distance measurements depending on the curvature prior assumed. 

fiducial model. This is automatically linked to the GUI and is accessible via the "F4C Extensions" 
drop-down menu in the GUI. 

• General growth function 

Growth constraints are computed exactly for any allowed w{z) and any curvature 17^. 

• Improved numerical derivative treatment 

The numerical derivative routine was modified from a simple double-sided derivative, to one which uses 
the complex-step algorithm [sot . with improved numerical stability, especially in extreme cases. 

• GUI button to include / exclude the prior information matrix 

A button has been added to the GUI to allow the user to either include or exclude the prior information 
matrix from the Fisher analysis. This makes it easy to compute and compare the impact of adding 
previous knowledge (e.g. Planck) on a given survey. 

E. User Extensions 

The general philosophy of Fisher4Cast was to make it as easy as possible to mould and extend to the needs 
of a general user. The power of Fisher4Cast does come at a price however and to modify the code requires 
the user to become familiar with the underlying structures used in the design and building of Fisher4Cast. 
The good news is that Fisher4Cast was designed to be as logical, elegant and general as possible and quick 
but limiting fixes were avoided where possible in favour of fiexibility. This means that once the user is 
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familiar with the underlying code structure, Fisher4Cast should be easy to extend in ways that were not 
even conceived of when it was written. 

All Fisher4Cast functions begin with the prefix FM while extensions have the EXT prefix (in ad- 
dition the figure plotting code is prefixed by FIG. To use Fishcr4Cast for a new application, code up 
your function (the X in the notation of Section , as a Matlab function (or functions if you have mul- 
tiple observablcs) and name them appropriately: FM_my .function. m. Then code the derivatives (i.e. 
dX/d0A, if you have them analytically, otherwise numerical derivatives will be used) and name appropri- 
ately, e.g. FM_my_derivs.m. Make sure all functions take in a data vector and base parameter vector 
(the fiducial model where the Fisher Matrix will be evaluated) as arguments and that they return a vec- 
tor, as in the header documentation of FM_function_l.m and FM_analytic_deriv_l.m. Then go to 
the input structure and replace the entry 'FM_analytic_deriv_l' in the field input. function_names 
with 'FM_my_function'. Then replace the field numderiv.f{l} — sprintf('FM_function_l'); with 
numderiv.fjl} = sprintf('FM_my .function'). 

As an example consider alternative parameterisations for dark energy. Currently, the Fishcr4Cast GUI is 
hard-coded for three cosmological observablcs {H,dA, and G), assuming the Chevallicr-Polarski-Lindcr (CPL) 
parameterisation of dark energy with parameters (z/;o,iOa), see Eq. ()lip . This is true of both the functions 
themselves, and the analytical derivatives included in the Fisher4Cast suite. The general framework of 
Fisher4Cast, however, means that one is not restricted to this parameterisation. As can be seen from 
Eq. p2|) . IHI) and ([22|) . the dark energy equation of state enters the cosmological observablcs through 
the evolution of the dark energy, via f{z), defined in Eq. (jl4p . Hence for any given 10(2) one only needs 
to specify the names of the functions (in the input structure) that will replace the current versions of 
FM_function_l.m {H{z)), FM_function_2.m (c?^(z)) and FM_function_3.m {G{z)). The same is true 
for the derivatives - either they can be coded analytically for the particular parameterisation of dark energy, 
or the derivatives will be evaluated numerically from the functions specified in the input structure. 

As a caveat, the GUI can only be used if the new parameterisation of dark energy still contains only two 
coefficients. If this is not the case, Fisher4Cast must be run from the command line. 

V. CONCLUSIONS 

The Fisher Matrix formalism is the standard forecasting method in cosmology and assuming one is 
familiar with it, allows rapid and widely understood results. In principle it is easy to learn and code for 
oneself. In practise there is a relevance threshold below which it simply is not worthwhile. 

Fisher4Cast facilitates both the acquisition of knowledge and the generation of production-quality Fisher 
Matrix forecasts. Written in Matlab, the Fisher4Cast Graphical User Interface (GUI) allows easy explo- 
ration of cosmological constraints and has both an interactive 'point-and-click' facility for automatic ellipse 
generation (when activated from the "F4C Extensions" drop-down menu) and an automatic M^t;X summary 
and results file generator, making direct inclusion of any output into research documents straightforward. It 
doesn't yet write the paper for you, but we are working on it. 

In this paper we have focused on illustrating novel uses of the Fisher4Cast suite by exploring the landscape 
of Fisher Matrix cosmology, as illustrated by Figures ([2])-®, as well as highlighting the effects of growth 
and curvature on Fisher Matrix forecasts of future cosmological surveys. These illustrate a limited set of 
applications of Fisher4Cast which we hope stimulates members of the community to use and extend the 
code. 

Research that has thus far used Fisher4Cast, either for forecasts or for producing plots, includes [5l|, [lH| 
and [11], as well as unpublished work for the ACT survey by RH. 
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FIG. 12; Varying redshift of H — examples of ellipses corresponding to particular values in the FoM landscape 
of Figure (O, in the particular case where the redshift bins for H are varied while keeping the bin for cIa fixed at 
2 = 0.6. The degeneracy direction for (Ia (black dashed line) and H (red solid line) are included to help see how 
their relative orientations contribute to the orientation and size of the combined ellipse. The degeneracy direction 
of H rotates anti-clockwise with the most rapid rotation experienced initially for the low redshift range of H. After 
z — 1 the rotation of H slows down and no longer is significant. As the H and dA degeneracy directions become 
more orthogonal so the resulting constraints improve, yielding a higher FoM. The ellipses are coloured corresponding 
to their FoM where red indicates the largest (around 0.08), and blue the smallest, FoM. 

The FoM landscape in Figure ([S]), and in particular the valleys separating high points in FoM space 
warrant further investigation. We argued previously that the variation in the FoM was an interplay between 
the widths of the ellipses (the ellipses for a single redshift measurement arc infinitely degenerate along the 
semi-major axis) and the orientation of the H and dA ellipses. This is illustrated in Figure where we 
show the two degenerate ellipses (red solid lines for H and black dashed lines for cIa) and the combined ellipse 
for a series where we keep the dA bin fixed at z = 0.6 while varying H in the redshift range 0.1 < z < 5. 
This moves us along a line in the landscape that intersects the peak at an H bin of z = 1.5. The combined 
ellipses have colours that correspond to the FoM in Figure ([5]), the highest value of the DETF FoM (^ 0.08) 
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coloured in red, and the lowest 0) coloured in dark blue. One can see a distinct rotation of the degeneracy 
direction for H as the redshift changes. This anti-clockwise direction rotation is most pronounced at low 
redshift from z = 0.1 to z = 1.0 after which it becomes more subtle and stops rotating completely at higher 
rcdshifts. 



APPENDIX B: FISHER DERIVATIVES FOR H{z),dA{z),G{z) 

We present the analytical Fisher derivatives of the Hubble parameter, H{z), given as Eq. and 
the angular diameter distance, (i^(z), given as Eq (jlSp with respect to the cosmological parameters 
(HQjfljm^kjWojWa), where WQ,Wa are the CPL dark energy parameters, and assume the forms f{z) and 
E{z) as given in Eqs. (fT4|) and (fT2|) respectively. In all cases the derivatives are taken in a general Friedmann- 
Lcmaitre-Robertson- Walker background without assuming flatness. 

1. The Hubble parameter 

As the Hubble constant, Hq only appears as a multiplicative term in H{z), the Fisher derivative of the 
Hubble parameter with respect to the Hubble constant Hq is simply 

ri ff 

Derivatives of the function S'{z) = H'^{z)/Hq = E'^{z) arc found in all derivatives of both H and dA and are 
worth defining separately: 



df{z) 

d£{z) 
dfl 



k 



{l + zY-f{z) 



{l + zY-f{z) 



— — — (i - ii„i - iifcj— — 

OWq OWq 

= 3(l-r!„-17fe)/(z)ln(l + z) 

OWa OWa 

= 3(1 - - nk)f{z) (^ln(l + z) - . (B2) 

(B3) 

For all the cosmological parameters we consider other than the Hubble parameter Hq, the derivatives with 
respect to the Hubble parameter can then be expressed as 

dHjz) Ho d£iz) . 
89- ^ 2E~de~' 0i ^ {^m,^k,wo,Wa) (B4) 

2. Angular Diameter Distance 

In a FLRW background, dA{z) is given by Eq (fT5|) . The function x(z) is defined by Eq. (fT6|) . It is again 
useful to define first order derivatives of x(z) with respect to the parameters, as these derivatives will occur 
many times in the Fisher derivatives of dA^z). 

dx(z) _ r 1 d'^i^'),, 



86, 



1 8(S'(z') 

Jo 2E^{-') 80- i^^'^k,u!o,Wa), (B5) 



20 

again where (a{z) = H'^{z)/Hq. The Hubble parameter appears only m the pre-factor of the angular diameter 
distance, and hence 



dHo 



1 c 1 

1 + y/Tk: 



sinh 



1 

i5o 



dA{z) 



(B6) 



The matter density Cl„i contributes solely to the xi^) term, and hence can be expressed using Eqs. (|B5p 
and (IB3I) as 
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with the derivatives of the co-moving distance given by 
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The curvature parameter is found both in the pre-factor and the sinh term of the angular diameter distance, 
hence 
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The Taylor series expansion of Eq (|B9[) is used Q,k 0, namely: 
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where 



X(Z,0) EE X{z) 



(Bll) 



are the functions (for example E{z), x{z)) assuming flatness. Using the definitions Eqs. (jBSP and ()B3p . the 
derivatives of the angular diameter distance with are expressed similarly for 9i € {wo,Wa) as: 
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3. Growth 



All derivatives related to the growth are computed numerically (in the routine FM_num_deriv.m) by 
solving Eq. (|22p and using either the complex step or central finite difference algorithm, depending on the 
choice of the user. 
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APPENDIX C: FISHER4CAST FIGURE CODE 
1. Code to produce Figure ([3]): plot of function values and derivatives 

% 

function FIG_f unction_derivative_plot (deriv_f lag, f unction_f lag) 
°L This function generates a plot of the derivatives of the specific 
°/o observables included in Fisher4Cast, H(data) , d_A(data) and G(data) . 
"/o It must be included in the directory in which Fisher4Cast is contained, 
"/o or that directory must be added to the Matlab path. 

°L The flags deriv_flag and f unction_f lag are set to 1 (0) if you do (don't) 
"/o want to plot the function or derivatives. If no input is given these are 
°L both set to 1 and you get plots of all functions and derivatives. 

"/o As a default example it calls the Seo_Eisenstein_2003 input structure, 

"/o but then generates a redshift vector from 0.1:10; 

°/» The colours for line plots must be specified as 1x3 RGB vectors, 

7o normalised to 1 (i.e. so each entry divided by 255). See the default 

"/o colours as an example . 

close all 

°L Flags to control what you want to plot, either derivatives only, of 
°L function only, or both 
if nargin == 

deriv_flag = 1; 

f unction_f lag = 1; 
elseif nargin == 1; 

f unction_f lag = 0; 

end 

I 

°L Specify the colours of the derivatives 
hcolour =[147 50 0]./255; 
gcolour = [240 201 81]./255; 
dacolour = [231 109 29]./255; 
colourmat = [hcolour 
gcolour 
dacolour] ; 
styles = {'-', ' — '-.'}; 

legendmat{l} = {MlnH/dH_0' ' dlnH/dln\Dmega_m' . . . 

'dlnH/d\Omega_k' 'dlnH/dw_0' 'dlnH/dw_a'}; 
legendmat{2} = {'dlnd_A/dH_0' ' dlnd_A/dln\Omega_m' . . . 

'dlnd_A/d\Omega_k' ' dlnd_A/dw_0 ' ' dlnd_A/dw_a' } ; 

legendmat{3} = {'dlnG/dH_0' ' dlnG/dln\Omega_m' . . . 

'dlnG/d\Omega_k' 'dlnG/dw_0' 'dlnG/dw_a'}; 
I 

°L Generate the Input data for the derivative plot 

input = Seo_Eisenstein_2003 ; Z initialise the input structure 

data = 0.1:0.1:10; 7o the redshift range we want to consider 

data = data( : ) ; 

input . growth_zn = ; 

input . growth_zn_f lag = 1; 

"/o make sure the Growth is normalised at data = 0; 

input . observable_index = [12 3]; °L Use all three observables 

input .num_observables = length(input . observable_index) ; 

°L Re-assign the redshift vectors in the input structure 

input . data{l} = data; 

input . data{2} = data; 
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input .data{3} = data; 

"/o Re-assign the errors vectors in the input structure 
input . error{l} = . 1 . *ones (1 , length(input . data{l}) ) ; 
input . error{2} = . 1 . *ones (1 , length(input . data{2}) ) ; 
input . error{3} = . 1 . *ones (1 , length(input . data{3}) ) ; 

"/o Use the analytical formula for H, d_A and 
"/onumerical derivatives for G 
input. numderiv.flag{l} = 0; 
input. numderiv.flag{2} = 0; 
input. numderiv.flag-CS} = 1; 

I 

"/o Run Fisher to get the parameter values and derivatives 
output = FM_run( input) ; 

close(l); "/o Close the figure of the Fisher Ellipse 

% 

t PLOT THE DERIVATIVES 

1 We will plot dlnX/dtheta_i = dX/Xdtheta_i 

I 

if deriv_flag == 1 

for i = 1 : input .num_observables 7. Loop over the observable functions 

X = input . observable_index(i) ; 
% Set the figure properties 
figure (x*100) 

axes( 'FontNcune', 'Times', 'FontAngle', 'italic',... 

'FontSize', 14 , 'XScale', 'log', 'XTickLabel' , {'0.1';'1';'10'}) 
hold on 
box on 

xlabelCRedshift' , 'FontName', 'Times',... 

'FontAngle', 'italic', 'FontSize', 16 ) 
ylabelC ['Fisher Derivatives for ', input . observable_names{x} ,.. . 

'(z)'], 'FontName', 'Times', 'FontAngle', 'italic', 'FontSize', 16) 

for j = 1:5 
if j==2 

7o The Omega_m derivative, this is actually dlnH/dlnOm 
semilogx(data, input .base_parameters (j ) . . . 

. *output . f unction_derivative{x}( : ,j) . /output . f unction_value{x} , . . . 
'LineStyle', styles{j}, 'LineWidth' , 2,... 
'Color', colourmat (x, : ) ) 

else 

plot(data, (output .function_derivative{x}( : , j) . . . 
. / output . f unction_value{x}) , . . . 

'LineStyle', styles{j}, 'LineWidth' , 2, 'Color', colourmat (x, :) ) 
end 7o end the check to see if we are plotting Omega_m derivatives 
end 7o end the loop over the parameters 

legend ( legendmat {x} , 'Location' , 'NorthWest') 
°L plot the legend for the function according to the observable 
end 7o end the loop over the observable functions 

end y„ end the if loop for plotting of derivs 



1 
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y. PLOT THE FUNCTIONS 

% 

if f unction_f lag ==1 

for i = 1 : input .num_observables 

X = input . observable_iiidex(i) ; 
% Set the figure properties 
figure (100*x +1) 

axes( 'FontNcune', 'Times', 'FontAngle', 'italic', ... 
'FontSize', 14 ,'XScale', 'log', 'XTickLabel' , { ' . 1 ' ; ' 1 ' ; ' 10 ' > ) 
hold on 
box on 

xlabelCRedshift' , 'FontName' , 'Times', 'FontAngle',... 

'italic', 'FontSize', 16 ) 
ylabel ( [input . observable_naines{x}- , '(z)'], 'FontName',... 

'Times', 'FontAngle', 'italic', 'FontSize', 16 ) 

7o Plot the data 

semilogx(data, output . funct ion_value{x} , 'LineWidth' , . . . 
2, 'Color', colourmat (x, : ) ) ; 

end 7o end the loop over observables 

end % end the if loop for function plotting 

2. Code to produce Figure (l6|): volume slice plots 

I 

"/oThis function takes a three dimensional matrix, fom_vol_out, and plots a 

7oSlice plot for this data. 

1 

"/of om_vol_out can be generated by calling 

y. 

y»>>f om_vol_out = FlG_generate_f om_volume_data(vol_res) 

y. 

yif no vol_res is passed a default value of 30 is assumed. 

y. 

y»Exajnple of using this function: 
I 

yo>>FlG_plot_slice_f om_volume(f om_vol_out) 

y. 

y.if no fom_vol_out is given then the code checks to see if there is a default 

y.mat file, def ault_f om_vol_out .mat , to load the data from, else the 

y.f unction FlG_generate_f om_volume_data is called with a set of default values. 
I 

function FlG_plot_slice_f om_volume (f om_vol_out) 

Zcheck if a matrix of the volume space, fom_vol_out, is passed to the function 
if nargin<l 

y.see if the default .mat file exists and load the data 
if exist ( ' def ault_f om_vol_out .mat ' ) 

load def ault_f om_vol_out ; 

fom_vol_out = def ault_f om_vol_out ; 

else 

'Lli not then generate the fom_vol data 

f om_vol_out = FlG_generate_f om_volume_data(30) ; 
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end 

end 

y„set the colormap 
colormap(jet) ; 

"/oselect the planes to intersect for the slice plot 

I 

"/othis is an additional subsection of code to make the selection of redshift 
"/oplanes to plot more generic and easy to manage for a range of 
"/of om_vol_out ' s produced. Please note you must still manually specify the 
y„redshift range 

[x_col y_col z_col] = size (f om_vol_out) ; 
x_redshif t_range = [0,5]; 
y_redshif t_range = [0,5]; 
z_redshif t_raiige = [0,5]; 

"/.calculate a relationship from column to redshift 

column_to_redshif t_ratio_x = x_col/x_redshif t_rELnge (end) ; 

column_to_redshif t_ratio_y = y_col/y_redshif t_range (end) ; 

column_to_redshif t_ratio_z = z_col/z_redshif t_range (end) ; 

y»specify the slices redshift to intersect the fom volume space 

x_redshif t_slice = [0.6667]; 

y_redshift_slice = [0.4167, 1.6667]; 

z_redshif t_slice = [2.5]; 

y.calculate the slices in column numbers 

xslice = column_to_redshif t_ratio_x . *x_redshif t_slice ; 

yslice = column_to_redshif t_ratio_y . *y_redshif t_slice ; 

zslice = column_to_redshif t_ratio_z . *z_redshif t_slice ; 

I 

yplot using slice 

s = slice(fom_vol_out, xslice, yslice, zslice) ; 

y.set the labels so they match the rsinge of 

y redshift as opposed to the column numbers 

set(gca, 'xtick' , [O:column_to_redshift_ratio_x:x_col] , . . . 

'xticklabel' , [0:x_redshift_range(end)] , 'ytick' , . . . 

[0 : column_to_redshif t_ratio_y : y_col] , 'yticklabel' , . . . 

[0:y_redshift_range(end)] , 'ztick' , [0 : column_to_redshif t_ratio_z : z_col] , . . . 
'zticklabel' , [0:z_redshift_range(end)] ) ; 

yset the X y and z labels 
ylabeK'H Redshift') ; 
xlabel('d_A Redshift'); 
zlabeK'G Redshift'); 



3. Code to produce Figure Fisher ellipses as a function of changing curvature prior 

I 

function FlG_vary_f om_curvature_prior 
global input plot_spec axis_spec 

y This function generates a plot of the Fisher ellipse as one changes the 
y prior value on curvature, eaid a corresponding plot of the Dark Energy 
y Task Force Figure of Merit (FoM) as a function of the prior, 
y See the User's Manual for definitions of the FoM. 
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7o This code must be included in the directory in which 
°L Fisher4Cast is contained, or that directory must be 
°L added to the Matlab path. 

"/o As a default example it calls the Seo_Eisenstein_2003 input structure, 

°L The Matlab 'colormap' command is used to generate the line colours, 

°L specific to each observable and the number of iterations is given by N. 

"/o NOTE That this code uses getfigdata.m by M.A. Hopcroft, 
"/o which is code from the Matlab File 

'/o Exchange (\protect\vrule widthOpt\protect\href {http : //www.mathworks . co.uk/matlabcentral/fileexchange/ 
°L It is included in this package. 

close all 

I 

7« Choose your colour schemes for the various combinations 

linecolor{l} = sort (colormap (copper) , 'descend'); 

linecolor{2} = colormap (autumn) ; 

linecolor{3} = colormap (copper) ; 

linecolor{4} = colormap (winter) ; 

linecolor{5} = colormap (summer) ; 

input . fill_f lag = 1; 

valinit = 35; °L the starting colour for the plots 
num_obs = [12]; 7. The vector of combinations you want: 
7. 1 = Hubble 

7o 2 = Angular Diameter distance 
7 3 = Growth Function 

7 4 = Hubble parameter + Angular Diameter distance + Growth 
7 5 = Hubble parameter + Angular Diameter distance 

line_width = 2; 

I 

7 Set the Input structure for the survey you will use 
input = Seo_Eisenstein_2003 ; 
input . data{3} = input . data{l}- ; 

input . error{3} = . 1 . *ones (1 , length(input . error{l}) ) ; 

input . observable_index = [12 3]; 7 We will use all observables 

input . fill_f lag = 0; 

input .numderiv. flag{3} = 1; 

I 

7 Set up the range you wish to consider 
start_prior= le6; 

range = le8; 7 No of orders of magnitude in the prior 
N = 20; 7 Number of points 
amp = (rajige) ~ (1/N) ; 

7 Initialise the priors 

prior_orig = input .prior_matrix; 7 this will be the default value 
input .prior_matrix(3,3) = start_prior; 
input .prior_matrix(2,2) = 0; 

7 Initialise the prior on the matter density to zero 

I 

7 Initialise the global FoM plot 
figure (3000) 

axes( 'FontName', 'Times', ' FontAngle ' , . . . 

'italic ' , ' FontSize ' , 14 , ' XScale ' , ' log ' , ' YScale ' , ' log ' ) ; 
hold on 
box on 
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xlabeK 'Prior (\Omega_k) ' , 'FontNaune' , 'Times',... 

'FontAngle', 'italic', 'FontSize', 16 ) 
ylabeK 'Figure of Merit', 'FontNaune', 'Times',... 

'FontAngle', 'italic', 'FontSize', 16 ) 



°L Call Fish.er4Cast in a loop 
for ni = num_obs(l) :num_obs(end) 

input .prior_matrix(3,3) = start_prior; 

if ni ==4 

°L Compute the Fisher Ellipse for combination of Hubble, d_A and G 

input . observable_index = [12 3]; 
elseif ni ==5 

°L Compute the Fisher Ellipse for combination of Hubble and d_A only 
input . observable_index = [12]; 

else 

input . observable_index = ni; °L use index value as specified 

end 

for i=l:N 

f igure(l) 
hold on 

input .prior_matrix(3,3) = input .prior_matrix (3,3) . /amp; 
7o modify the Prior 

output = FM_run(input) ; °L Call Fisher4Cast 
val(i) = input .prior_matrix (3,3) ; 
7„ save the value of the prior for plotting 
outv(i,:) = output. fom; 

°L Save the full FoM vector 
out(i) = outv(i,l); °L Save the DETF FoM 

h = getf igdata(l) ; 

'/o call getfigdata.m to rip off the ellipse 

x{i} = h{l}.x; 

y{i} = h{l}.y; 

close (1) 7o close the figure 

end 

I 

°L Initialise the plotting figures with the axis specs etc 
figure (1000+ni) 

axes( 'FontName', 'Times', 'FontAngle', 'italic', 'FontSize', 14 ) 
hold on 
box on 

xlabel ( 'w_0 ' , 'FontName', 'Times', 'FontAngle',... 

'italic', 'FontSize', 16 ) 
ylabel ( 'w_a' , 'FontName', 'Times', 'FontAngle',... 

'italic', 'FontSize', 16 ) 
axis([-3 1 -10 10 ]) 
count = ; 

I 

°L Plot the resulting ellipses 
for i = 1:N 

f igure(1000+ni) 

hold on 

°L Use increasing or decreasing colour to get a gradient 
if i < ceil(N/2) 

count = count +1; 

else 

count = count -1; 
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end 

cval = valinit+ 2*count; 

plot(x{i}, y{i}, 'Color', liiiecolor{iii}(cval, : ) , . . . 
'LineWidth' ,line_width) 

end 

figure (3000) 
hold on 

plot (val , out , 'Color', linecolor{ni}-(20, : ) , . . . 
'LineWidth', line_width) 

end 

I 



APPENDIX D: QUICK START GUIDE TO FISHER4CAST 

1. Hardware and software requirements 

This software is written to be run in Matlab (Linux, Windows and under Mac OSX, although this has not 
been extensively tested). The user needs Matlab installed (Tested on Version 7) to be able to run this code. 
Free disk space of approximately 2MB and the minimum recommended processor and memory specifications 
required by the Matlab version you are using is suggested. 

2. Downloading Fisher4Cast 

Currently the code is available for download at one of the following websites [H, [s^l. Save this .zip file 
into the directory you want to run the Fisher4Cast suite from. 

3. Getting started 

The code can be run from the command line or the Graphical User Interface (GUI). Wc describe the 
command line below, and mention how to get the GUI started. More information on the GUI can be found 
in the Users' Manual [ssj . 

4. The Graphical User Interface 

• Running the GUI 

The GUI can be started from the Matlab editor. The file FM_GUI.m must be opened from the 
directory, and once the file is opened (click on the file icon from within the Command-line interface to 
open it with an editor) press F5 to run the code. This will open up the GUI screen. 

You can also launch the GUI from the command line by typing: 
»FM_GUI 

This then functions in the same way as using FM_run in the command line (as explained in the following 
section). 

For more information on the technicalities of the GUI, see the full manual. 

• GUI Screenshots 

We include some screenshots of the Graphical User Interface. 
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FIG. 13: Plotting multiple ellipses on one axis - using the 'Hold on' multiple error ellipses can be overlaid on 
one axis. The 'Area Fill' command allows you to choose the colours for the error ellipses. Also shown is the 'Running' 
window which indicates the code is running to calculate the Fisher ellipses. 




FIG. 14: Different background images and colour schemes - the background images and colour schemes (skins) 
allow for a fully customisable Graphical User Interface. 




FIG. 15: Various Figures of Merit can be plotted - the drop-down list allows for a choice between various 
Figure of Merit options. 



5. The Command Line 



• Running the code 
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Open your version of Matlab and change the working directory to be the same as where you saved 
Fisher4Cast in. To run the code from the command hnc with one of the standard test input structures 
suppHed, type: 



»output = FM_run(Cooray_et_al_2004) 



This will call the code using the pre-supplied test input data (Cooray_et_aL2004) and then generate 
an error ellipse plot for the parameters and observables supplied in the chosen input. All the relevant 
generated output is written to the output structure. You can sec the range of outputs to access by 
typing: 



>>output 



and then examine each output individually by specifying it exactly. For example: 



>>output . marginalised_matrix 



will access the marginalised Fisher Matrix from the output structure. 

You can use the supplied input files as a template for generating new input files with your own 
customised parameters and values. All fields shown in the example structures must be filled in any 
user-defined structure. 

The code can also be run from the Matlab editor. Once the code is opened (open it from inside the 
Matlab window), you can press F5 to run the code. Note that if the code is run from the Editor it 
will call the default input structure, which is the Cooray_et_al_2004.m file. This is an example file 



containing input data from the paper by Cooray et al. |54l | . This output can be directly compared to 
that of Figure 1 of that paper. If your output compares correctly you have a working installation of 
the code. Another input available is Seo_Eisenstein_2003.m [381 . 



